Graviton production in non-inflationary cosmology 
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We discuss the creation of massless particles in a Universe, which transits from a radiation-dominated 
era to any other expansion law. We calculate in detail the generation of gravitons during the tran- 
sition to a matter-dominated era. We show that the resulting gravitons generated in the standard 
radiation/matter transition are negligible. We use our result to constrain one or more previous 
matter-dominated era, or any other expansion law, which may have taken place in the early Uni- 
verse. We also derive a general formula for the modification of a generic initial graviton spectrum 
by an early matter dominated era. 
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I. INTRODUCTION 

One of the most interesting aspect of inflation is that 
it leads to the generation of a scale-invariant spectrum 
of scalar perturbations [l[ and of gravitational waves Q 
(see also 0, H). The origin of these perturbations is 
the quantum generation of field correlations in a time- 
dependent background (scalar field modes for the scalar 
perturbations and gravitons for the tensor perturba- 
tions). Since this generation takes place mainly on super- 
horizon scales, it is not correct to talk of 'particles'. How- 
ever, long after inflation, when the perturbations re-enter 
the horizon, the particle concept, e.g. for gravitons be- 
comes meaningful and we can calculate, e.g. the energy 
density of the gravitons which have been generated dur- 
ing inflation. 

It is natural to ask whether particle production takes 
place also in an ordinary expanding but non-inflationary 
Friedmann Universe. The answer is that particle produc- 
tion (or more gcncrically the quantum generation of field 
correlations) can indeed take place after inflation, but 
if there is no inflationary phase to start with, the initial 
vacuum state is in general not known, and the production 
rate cannot be computed. Examples where particle cre- 
ation taking place after inflation (or after pre-big bang) 
modifies the final spectrum are given in Refs. [H, [<!, 0, HI • 
Especially in Ref. [7[ it has been studied how inflationary 
perturbations are modified if the subsequent expansion is 
not standard radiation but some other expansion law. 

In general, the vacuum state, and hence the particle 
concept, is well defined only if the spacetime is static or 
very slowly varying || . Let us consider a mode of fixed 
(comoving) frequency k in a Friedmann universe. The 
above condition then corresponds to k/H 3> 1, where 
H = aH is the comoving Hubble scale. In this sense, the 
scale (wavelength) under consideration must be "inside 
the horizon" . However, the production of a particle with 
a given energy k can only take place if the energy scale 
of expansion is larger or of the order of the energy of 
the associated mode, i.e. k <, Ti . Therefore, having a 
well defined initial vacuum state, and subsequent particle 



creation, usually requires a decreasing comoving Hubble 
rate. This is verified only during inflation or during a 
collapsing Friedmann Universe, like in the pre-big bang 
or in bouncing models. 

However, one important exception to this general rule 
exists, and it is the subject of the present paper: in 
a radiation-dominated Friedmann background, massless 
perturbations do not couple to the expansion of the Uni- 
verse, and evolve like in ordinary Minkowski space. This 
has already been realized and studied to some extent in 
Ref. In a radiation-dominated Universe we there- 

fore can provide vacuum initial conditions for all modes 
of a massless field, including super-horizon modes. Thus, 
when the expansion law changes, e.g. from radiation- 
to matter-dominated, the massless modes couple to the 
expansion of the Universe, and those with k/H < 1 are 
amplified. 

In this paper we study this phenomenon in two situ- 
ations of interest. In the first, we investigate graviton 
production during the ordinary radiation/matter transi- 
tion at redshift z cq ~ 3500. We determine the amplitude 
and the spectrum of the generated gravity wave back- 
ground, and we show that the spectrum is flat and the 
amplitude is negligibly small. In the second case, we in- 
vestigate the production of gravitons during an arbitrary 
matter-dominated phase, which could take place in the 
early Universe, e.g., if a (very weakly interacting) particle 
becomes massive, and succeeds to dominate the Universe 
for a period of time before it decays into radiation. We 
derive a general formula for the gravity wave spectrum 
generated by any number of such intermediate periods 
of matter domination. We also determine the gravity 
waves produced by a transition into an arbitrary other 
expansion law. Finally, we discuss the modifications of 
our results which occur when the initial state is not the 
vacuum but some arbitrary state which may already con- 
tain particles. In this work we concentrate on graviton 
production, but all our results are equally applicable to 
other massless particles. 

The reminder of this paper is organized as follows. In 
the next section we present the setup and the basic for- 
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mulae used in our work. In Section III, we calculate the 
gravitational wave production during the standard radi- 
ation/matter transition. We also give the results for the 
transition from radiation to some generic expansion law. 
In Section ITVl we consider the effect of one or several ad- 
ditional transitions in the early Universe and we derive 
results for general non- vacuum initial conditions. In Sec- 
tion [V] we derive the consequences of our results and we 
draw some conclusions. 

Notation: We work in a spatially flat Friedmann Uni- 
verse, and we denote conformal time by t, so that 



ds 2 = a 2 (t) (-dt 2 



Sijdx l dx J 



An over-dot denotes the derivative with respect to the 
conformal time. We use natural units c = % = 1, except 
for Newton's constant G, which is related to the reduced 
Planck mass by 8ttG = m~ 2 . We normalize the scale 
factor, so that ao = a(to) = 1 at the present time. 



II. GRAVITON CREATION IN COSMOLOGY 

We now consider tensor perturbations of the Fried- 
mann metric, namely 

ds 2 = a 2 {t) [-dt 2 + (S i:j + 2h ij )dx i dx j ] , 

where hij is a transverse and traceless tensor. In Fourier 
space we have 

MM) = MMh ( + } (k) + MM)i _) (k) , (i) 

where (k) denote the positive and negative helicity 

polarization tensors, and fe*e^ =0. In a perfect fluid 
background, i.e. if there are no anisotropic stresses, both 
amplitudes satisfy the same wave equation, 



Eq. ([4} is the action of a canonical scalar field with time- 
dependent effective squared mass m 2 (t) = —a/ a. If the 
expansion of the Universe is slow enough (compared to 
the frequency of the mode under consideration) , then the 
effective mass is negligible, and the theory describes a 
massless scalar field in Minkowski space, which can be 
quantized according to the usual procedure: we first pro- 
mote the field to an operator 



(x,f) 



d 3 k 



e lk - x Xk(<)&k + e- ik -*xUt)bt] > (6) 



(2tt) 3 

then we impose the commutation rules 



[& k ,6 k ,] = (2^ 3 (k-k') 



0. 



The field equations derived from the action lead to 
the mode equation 



Xk 



k 2 -- 
a 



Xk = 0. 



(7) 



Within linearized gravity we can therefore quantize the 
metric fluctuations, provided the Universe expands adi- 
abatically, by making use of the above rescaling of the 
amplitude h. 

In particular, we now assume that the Universe is 
initially radiation-dominated, so that a(t) = t, and 
d/a = 0, and <j> represents exactly a massless scalar field 
in Minkowski space. We consider the vacuum initial con- 
ditions for the modes Xk(i) as given by 



Xk(t) = 



1 



-ikt 



(8) 



More general initial conditions will be considered at the 
end of Sec. IV. The field normalization is determined by 
the Klein- Gordon norm 



Oh = h + 2-h + k 2 h = Q , 
a 



(2) 



iX* d x = i{x*dox - xdoX*) = 1 



(9) 



where h = h±. This equation of motion is obtained when 
expanding the gravitational action in a Friedman uni- 
verse to second order in h, 

in 2 f 

S + SS = — ^ / d*xy/-(g + 6g)(R + SR) . (3) 

A brief calculation shows that the lowest (second order) 
contribution to 6S can be written in Minkowski-space 
canonical form, 

5S^=-\ Jd^x (d, 

if we rescale h as 

fc(x,t) 



09^0- -<f>' 

a 



J d A xC, (4) 



— — < 

V2 m p a(t) 



(x,i) 



(5) 



Then, the field operator (f> and its canonically conjugate 
momentum, II = d£/d(do(j)) satisfy the canonical com- 
mutation relations. The operators 6k define the in vac- 
uum by 6k|0in) = Vk. In the following, this is the initial 
vacuum, void of particles by construction. 

Suppose now that at a time t = t\, the Universe 
changes abruptly from radiation-dominated to another 
expansion law. Then, the effective squared mass no 
longer vanishes, and the initial modes, with kt\ < 1 are 
amplified. Continuity requires that \ an d X match at 
t = ti, and these conditions determine the Bogoliubov 
coefficients, which relate the new modes Xout and opera- 
tors 5k, out to the old ones, x; n and Sk,i n by @ 



Xout(i) = a-Xin(i) +/3x* n (t) , 



(10) 



(11) 
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With these relations, we can compute the number density 
of the particles [l| created at the transition Q 



iV(Mi) HOi„l&LuA,out|o in ) 



(12) 



Thus, the energy density p = (T o) can be written as 



1 



Ph 



d 3 k 



kN(k,h) 



1 



(2tt) 3 v ' ' 2ir 2 a 4 
which implies the usual formula [b3] 



dkk A 



dph 



d\og k 



, (13) 



(14) 



where we have multiplied Eq. (|13p by a factor 2 to take 
into account both polarizations. Note also that k denotes 
comoving momenta/energy so that we had to divide by 
a 4 to arrive at the physical energy density. The second 
quantity of interest is the power spectrum Ph(k,t), de- 
fined by 



r flu 

J y^( fc ^) = (°in|^(i>x) 2 |O il 

Using Eqs. ([5]) and ©, we obtain 

k 2 \3\ 2 



P h (k,t) 



{2nfm 2 p a 2 ' 



(15) 



(16) 



where, again, we have multiplied by 2 to account for both 
polarizations. 

Note that for all this it is not important that we con- 
sider a spin 2 graviton. The exactly same mode equa- 
tion is obtained for a scalar field and also for a fermion 
field. In the latter case, the commutation relations have 
to be replaced by the corresponding anti-commutation 
relations. 



III. FROM RADIATION TO MATTER ERA 

Before discussing a transition from the radiation- 
dominated era to the matter era, let us consider the tran- 
sition from radiation to some generic power law expan- 
sion phase, a oc t q with q ^ 1 at some time t\. In the new 
era a/a = q(q — l)/t 2 ^ 0. Note that only if q > 1 (or 
if q < which corresponds to inflation or contraction), 
so that m 2 (t) — —q(q — l)/t 2 is negative, we will have 
significant particle production. Since the expansion law 
is related to the equation of state parameter w = P/ p 



9 : 



1 + 3w 



(17) 



this requires w < 1/3. 

Let us start in the vacuum during the radiation era, 
then x is given entirely by the negative frequency modes, 
Eq. ([5]) . This means that we consider the situation where 



there are no significant gravity waves present from an ear- 
lier inflationary epoch. Here, we really want to study the 
production due solely to the radiation/matter transition. 
The general solution of the mode equation ([7]) in the new 
era are the spherical Hankcl functions [ll| of order q — 1 , 



Xk(*) = 



a i 



zh 
2k q 



(2) 



Pi 



(18) 



where z — kt. Note that inside the horizon, i.e. for 
z 1, zh ' x {z) oc exp(— iz) corresponds to the negative 

frequency modes while zh^l^z) oc exp(iz) corresponds 
to positive frequency modes. We match x and x at t = t\ 
to the radiation-dominated vacuum solution ((SJ) . A brief 
calculation yields the coefficients {z\ = kti) 



Oil 

0i 



-r 



(zz 1 +g)^ 1) 1 (z 1 )-z 1 / l «(z 1 )] , (19) 
(izi + q)hf\(z 1 )-z 1 h^(z 1 )] . (20) 



This instantaneous matching condition is good enough 
for frequencies for which the transition is rapid, i.e. z\ ^ 
1. In fact, for frequencies with z\ > 1, the transition is 
adiabatic and no particle creation will take place. This 
can also been seen when considering the limits of the 
above result for large z\. Then ai -> 1 and (3\ — > 0, but 
strictly speaking the above approximations are not valid 
in this regime where no particle creation takes place. We 
therefore concentrate on Z\ <C 1. 

Let us now study the specific case of the radiation- 
matter transition, i.e. q = 2. Then we have to consider 
spherical Hankel functions of order 1 and the solution is 
given by 



Xk(t) 



ai z - 



Pi z + i 



2k 



The matching at t\ now yields 
1 



1 



i 

Zl 



2z{ 



2 ! 



pi = 



2k 



2z 2 



-1iz\ 



(21) 



(22) 



We want to evaluate the quantum field <j) at late time, 
when z»l and the mode k under consideration is sub- 
horizon. Then, the solution (f2Tj) is again the Minkowski 
solution, 



Xk{t) ^ ^=e 



Pi 



2k 



2k 



(23) 



The number of gravitons generated during the matter 
era (before z ^> 1) is, see Eq. (TI"2"|) 



N(k,t) = |ft| s 



(24) 



The graviton power spectrum is given by Eq. (|16p . and 
the energy density by Eq. (pi 
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Using that p lax ia 4 



Prad(*l)af 



\nip-L\a\ and 



|/?i| 2 = z 1 4 /4, we obtain 

dfi h (t) _ 2f* rad k 4 



dlogfc 



3tt 2 m?W^ 



= n 



rati 



ffoft 



36 x 30 V m 



m 2 



Ti 



^rad /" Hi 

m p 



6tt 2 



(25) 



For the second equal sign we have used that TL\ = X/t\ = 
a\Hi, which is strictly true only in the radiation era, in 
the matter era we have TI = 2/t and at the transition a 
value between 1 and 2 would probably be more accurate. 
But within our approximation of an instant transition, we 
do not bother about such factors. For the last equal sign 
we used H\ = /0 ra d/(6m^) with p rad = ^f^T 4 where 
9eS = Nb + -§Nf is the effective number of degrees of 
freedom. 



For a generic transition we obtain |/3i 



12 _ -2q 



that 



dn h {k) 

d\og k 



rad 



4-2 9 



/4 so 



(26) 



This spectrum is blue (i.e. growing with k) if q < 2 and 
red otherwise. 

As z\ < 1 in the regime of validity of our formula, we 
need a red spectrum i.e. q > 2 to enhance the gravi- 
tational wave energy density with respect to the result 
from the radiation — > matter transition. According to 
Eq. (JT7J), this requires — l/3<w<0, a slightly negative 
pressure, but still non-inflationary expansion. 

In the standard radiation to matter transition when 
three species of left handed neutrinos and the photon are 
the only relativistic degrees of freedom, we have <? e ff = 
29/4. For this transition 



eq 



0.85eV : 



0.35 x 10~ 27 r 



'■p , 



hence the result (|25|) is completely negligible. 



IV. MORE THAN ONE RADIATION-MATTER 
TRANSITION 

We now consider an early matter dominated era. At 
some high temperature T\ T cq , corresponding to a co- 
moving time t%, a massive particle may start to dominate 
the Universe and render it matter-dominated. At some 
later time t 2 , corresponding to temperature T 2 , this mas- 
sive particle decays and the Universe becomes radiation- 
dominated again, until the usual radiation-matter tran- 
sition, which takes place at T eq = T 3 . We want to deter- 
mine the gravitational wave spectrum and the spectral 
density parameter dQf l /dlog(k) as functions of T\ and 

Let us first again start with the vacuum state in the 
radiation eta before t±. When, we just obtain the results 
for the Bogoliubov coefficients a\ and f3\ after the 



first transition. To evaluate the matching conditions at 
the second transition, matter to radiation, we set 



a 2 _ e -ikt , 02 c ikt 



2k 



t>t 2 . 



Matching \ and X & t ^2 we can relate the new coefficients 
a 2 and j8a to a\ and /3%. A brief calculation gives 



ai 1 



f 



Q 2 



A = A 1 + 



0i 



z 2 2z\ 



&if(z 2 ) + P\g{z 2 ) , 
i 1 



z 2 2zo 



Pif(z 2 ) + uig{z 2 ) 



2z\ 



a l - 2 iz 2 

2zf 



or in matrix notation 



« 2 



M(z) 
M-\z) 



= M(z 2 ) 



Oil 

01 



with 



/(*) 9(z) 

m i(z) 

m -g(z) 

-m m 



(27) 
(28) 
(29) 



(30) 

(31) 
(32) 



The fact that M € Sl(2, C), i.e., \ f(z)\ 2 ~ \g{z)\ 2 = I en- 
sures that the normalization condition §§§ which trans- 



lates to the condition 



= I for the Bogolioubov 



coefficients of a free field, is maintained at the transition. 
Finally, the matching at the usual radiation-matter tran- 
sition yields 



03 



a 2 1 



i 1 

z 3 2zl 



h[ '-J 3 -2zl)-2zl e 



02 p 2iz 3 

2zf 

"2 -2iz: 



Oi 3 
03 

a 3 

03 



a 2 
'2 



M- 1 (z 3 )M(z 2 )M- 1 (zi) 



(33) 
(34) 
(35) 



To obtain the power spectrum and energy density in 
this case, we simply have to replace |/3| 2 in Eqs. (fT6|) and 



JE} by |/? 3 | 2 . In Fig. Q] we plot \/3 3 \ 2 as a function of 
z 3 for different choices of t 2 . The instantaneous transi- 
tion approximation breaks down for z 3 > 1, hence only 
the left side of the vertical line is physical. For the right 
side one would have to solve the mode equation numeri- 
cally, but since we know that particle production is sup- 
pressed for these frequencies, we do not consider them. 
We concentrate on 23 < 1. For these wave numbers, also 
z\ < z 2 < z 3 < 1. 
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k 

FIG. 1: The Bogoliubov coefficient \f3'i(k)\ 2 for ti = 1, and 
t$ — 2000, with various values of ti- Only the k- values left of 
the vertical dashed line satisfy ktz < 1. They show clearly a 
k~ 4 slope and the amplitude is well approximated by (|38[) . 



back to radiation at time t 2 ni 1 < n < N , with the result 



\0: 



2AT+1 1 



T$ ■ ■ ■ T 2 n-iT c < 



T 2 ---T- 



2N 



(41) 



Hence, each return to the radiation-dominated era at 
some intermediate temperature T 2n leads to a suppres- 
sion factor (T 2n+ i/T 2n ) 2 , where T 2n +i denotes the tem- 
perature at the start of the next matter era. 

On large scales, kt cq < 1, the energy density spectrum 
is fiat. The best constraints on an intermediate radiation- 
dominated era therefore come from the largest scales, i.e. 
from observations of the cosmic microwave background 
(CMB) as we shall discuss in the next section. 

We now briefly consider the case when the initial con- 
ditions differ from the vacuum case, Eq. (JSJ). We assume 
an arbitrary initial state of the field given by 



Xk(t ) = 4Le- fc * + -S=e ifct 



2k 



2k 



(42) 



together with the normalization condition which ensures 
that the field is canonically normalized, |ao| 2 — \Po\ 2 = 1- 
The same calculations as above now yield 



This allows the following approximations, 



«3 - /?3 



1 

2z} 

iz 2 

2z 2 
'3z 2 z 3 



ai ~ Pi ~ 2 , Z\ <C 1 



a 2 ~ -fj 2 — — , z 2 < 1 



, z 3 < 1 . 



(36) 
(37) 
(38) 



To obtain the results ([57)1 and (|3"5|) we have to expand the 
exact expression (|27|) to fourth order and (|33|) to second 
order, but we consider only the largest term in the result 
given above, using also Zi < z 2 < z 3 . Therefore, in the 
approximate expression (f3"Tj) . where we have neglected 
a term proportional to Zi/z 2 , one no longer sees that 
a 2 — > 1 and /3 2 — > when t 2 — > ti and hence z 2 —* Zi- In 
this case there is no intermediate matter-dominated era 
and therefore no particle creation, hence /3 2 = 0. This 
can be seen from the exact expression given in Eq. (|2"T|) . 
Within these approximations, Eqs. (fl7q|) and ([14")) lead 

to 



Ph(k) 
dlogk 



1 



to 



(27r) 3 \m p atft 3 k 

flrad fe 4 |/3 3 | 2 
7T 2 afPrad(il) 

geffPl) / 7j_ 

18 x 45 \ m„ 



kU < 1 



(39) 



ad 



T 

cq 

i7 



kt ca < 1 



(40) 



This result can be generalized to several, say N, interme- 
diate radiation — > matter transitions at times t 2n -i an d 



0'3 



Af- 1 (z 3 )M(z 2 )M" 1 (z 1 ) f ^ 



(43) 



where M{z) is the matrix giving the transition from mat- 
ter to radiation defined in Eq. (f3Tj) . 

Expanding this in Zi , z 2 and Z3 , using Zi < z 2 < z 3 < 1 
one finds that to lowest non-vanishing order, the final 
result for fl 3 depends only \ctQ + (3q\- However, if the 
phase of «o and (3q are nearly opposite, i.e., ao — — flo > 
and if |ao| and therefore also \{3o\ are much larger than 1, 
a correction proportional to | cto — /?o | becomes important. 
More precisely, the last of Eqs. (|3"6"|) now is replaced by 



2z 2 
3ziz 3 



— {a a + f3 ) + 2i(a - Po) 
zi 



(44) 



If ao = 1 and (3q = 0, the second term can be neglected 
with respect to the first one and we reproduce the pre- 
vious result (|38p. As we see from this equation, a large 
phase difference between ao and j3o changes not only the 
amplitude but also the slope of the spectrum. Of course 
in concrete examples, like for a previous inflationary pe- 
riod, see Ref. 0], the coefficients ao and @o also depend 
on the wavenumber. 

In Fig. [5] we show the dependence of \P 3 \ 2 on |ao| for 
different values of the relative phase between ao and (3q 
(top panel) and as a function of the relative phase for dif- 
ferent values of |ao|. The difference of \(3 3 \ 2 between the 
case where ao and /3o are perfectly in phase and of op- 
posite phase is of the order of l/zi, if |ao| is significantly 
larger than 1. This is already evident from Eq. (|44[) . In 
the Fig. [2] we have chosen Z\ =0.1, a unrealistically high 
value, in order to have a better visibility of the phase 
dependence which then changes \P 3 \ 2 only by one order 
of magnitude. 
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FIG. 2: In the top panel we show |/?3| 2 as a function of «o for 
ao and f3o in phase (top, solid, black line), of opposite phase 
(lowest, dashed, red line) and with a phase difference of 0.8-7T 
(middle, dot-dashed, blue line). 

In the bottom panel we show |/?3| 2 as a function of the phase 
difference ip for |ao| = 80 (top line) |ao| = 40 (middle line) 
and |ao| = 10 (lowest line). We have chosen z\ = 0.1 in 
this plot and the overall vertical normalization is in units of 

(2z 2 /(ziz 3 )) 2 . 



In conclusion, in the case of a non- vacuum initial state, 
|ao| significantly larger than 1, graviton production is en- 
hanced typically by a factor of order |ao| 2 ~ |/3o| 2 which 
is the number of initial particles. Hence in addition to 
the spontaneous creation we now also have induced par- 
ticle creation which is proportional to the initial particle 
number and much larger than the spontaneous creation 
if the particle number is large. An interesting point is 
that the phase shift between ao and (3q can significantly 
affect the final spectrum. 



V. DISCUSSION AND CONCLUSIONS 

The fact that the observed CMB anisotropics are of the 
order of 10 -5 yields a strong limit on gravitational waves 
with wave numbers of the order of the present Hubble 
scale, see e.g. [i~2j . 



d log k 



< 10" 



(45) 



k=H 



With f2 ra d ~ 10 , and Eq. (|40|). this implies the limit 



II) ^a<io- 5 . 



in. 



To 



(46) 



We know that during nucleosynthesis the Universe was 
radiation-dominated, hence T 2 > O.lMeV. With T cq ~ 
lcV, the above inequality reduces to T\ < m p for the 
value T 2 ~ O.lMeV, and it is even less stringent for 
higher values of T 2 . Hence even though the production of 
gravitons during an intermediate matter era is of princi- 
pal interest, we cannot derive stringent limits on T\ and 
T 2 . On the other hand, for values of T\ and T 2 close to 
the maximal respectively minimal value, T\ ~ m p and 
T 2 ~ O.lMeV, these gravitons would leave a detectable 
signature in the cosmic microwave background. 

We can, however, use this effect to limit any intermedi- 
ate era with q > 2, i.e. — 1/3 < w = P/p < 0. According 
to Eq. (|2"6"1) , in the general case, the particle number is of 



the order of \/3i\ 2 ~ z t 2q ', so that 



dlog k 



^rad(fc*l) 



4-2(j 



T 

1 cq 
~T~2 



z x < 1 . 



As above, T 2 denotes the temperature at which the Uni- 
verse returns to the radiation-dominated state, hence 
nucleosynthesis requires T 2 > O.lMeV. In this case, if 
q > 2, the spectrum becomes red and, at k ~ Hq, the 
limit can become quite interesting. Hence, graviton pro- 
duction can significantly limit a (non-inflationary) phase 
with negative pressure in the early universe. For k = H a , 
using t\ = Tii 1 = 1/ '(ai-ffi), we obtain 



(floti)" 



geffpl) 
10 5 



1/2 



where T ~ 0.2 x 10 _4 eV is the present temperature of 
the Universe. This can be a significant factor for large 
values of T%. Inserting this expression in Eq. (I47p . the 
limit (gSJ yields 



9os{Ti) 
10 5 



9-2 



6(9-2) 



m, 



,10" 



-10 



For example, for w 
obtain 

Ti 
T 2 



•1/21, (q 



m„ 



10- 



T 

± cq 

~T 2 ~ 

(48) 

7/3), we already 



(49) 
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E.g. for T 2 = lMeV this implies T x < 10 8 GeV. For 
smaller values of w the limit becomes more stringent. 

In this paper we have shown that there is cosmo- 
logical particle production of massless modes whenever 
the expansion law is not radiation-dominated so that 
a/a =/= 0. This term acts like a time-dependent mass 
and leads to the production of modes with comoving 
energy k 2 < \d/a\ ~ Ti. 2 , hence with physical energy 
uj — k/a^H. One readily sees that particle production is 
significant only if q < 1, i.e. the squared mass ~a/a < 0. 
When starting from a radiation-dominated Universe, we 
do have a well defined initial vacuum state also for the 
super-horizon modes which can be amplified by a tran- 
sition to another expansion law. We have also given the 
expressions for the produced particle number in the case 
of a generic, non- vacuum initial state, Eq. (p4"4"l) . This 
can be applied to an arbitrary inflationary, pre-big bang 
or bouncing model, which may already contain gravitons 
before the first radiation era. 

We have explicitly calculated the production of gravi- 
tons for a vacuum initial state and have arrived at the 
following main conclusions: 



i) The gravitons produced after the standard 
radiation-matter transition are completely negligi- 
ble. 

ii) If we introduce a matter-dominated era in the early 
universe, this leads to a flat energy spectrum of 
gravity waves which, in the most optimistic case, 
can be sufficient to contribute to the CMB tensor 
anisotropies in an observable way. 

iii) A phase of expansion with — l/3<w<0 leads to 
a red spectrum of gravitons. Such a phase in the 
early universe is severely constrained mainly by the 
amplitude of CMB anisotropies. 

iv) A graviton spectrum present at the beginning of 
the radiation era can become significantly amplified 
and modified by intermediate not-standard evolu- 
tion of the universe. 
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